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Abstract 

We present a general procedure for obtaining the coefficients of the scalar bubble and triangle 
integral functions of one- loop amplitudes. Coefficients are extracted by considering two-particle and 
triple unitarity cuts of the corresponding bubble and triangle integral functions. After choosing 
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unconstrained parameters of the cut loop momentum approach infinity. In this way we can produce 
compact forms for scalar integral coefficients. Applications of this method are presented for both 
QCD and electroweak processes, including an alternative form for the recently computed three- 
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of this extraction procedure allows for a very straightforward automation of the procedure. 
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I. INTRODUCTION 



Maximising the discovery potential of future colliders such as CERN's Large Hadron 
Collider (LHC) will rely upon a detailed understanding of Standard Model processes. Dis- 
tinguishing signals of new physics from background processes requires precise theoretical 
calculations. These background processes need to be known to at least a next-to-leading 
order (NLO) level. This in turn entails the need for computation of one-loop amplitudes. 
Whilst much progress has been made in calculating such processes, the feasibility of produc- 
ing these needed higher multiplicity amplitudes, such as one-loop processes with one or more 
vector bosons (W's, Z's and photons) along with multiple jets, strains standard Feynman 
diagram techniques. 

Direct calculations using Feynman diagrams are generally inefficient; the large number 
of terms and diagrams involved has by necessity demanded (semi) numerical approaches be 
taken when dealing with higher multiplicity amplitudes. Much progress has been made in 
this way, numerical evaluations of processes with up to six partons have been performed [1, 
2, 3, 4, 5]. On assembling complete amplitudes from Feynman diagrams it is commonly 
found that large cancellations take place between the various terms. The remaining result is 
then far more compact than would naively be expected from the complexity of the original 
Feynman diagrams. The greater simplicity of these final forms has spurred the development 
of alternative more direct and efficient techniques for calculating these processes. 

The elegant and efficient approach of recursion relations has long been a staple part of the 
tree level calculational approach [6, 7]. Recent progress, inspired by developments in twistor 
string theory [8, 9], builds upon the idea of recursion relations, but centred around the use of 
gauge-independent or on-shell intermediate quantities and hence negating a potential source 
of large cancellations between terms. Britto, Cachazo and Feng [10] initially wrote down a set 
of tree level recursion relations utilising on-shell amplitudes with complex values of external 
momenta. Then, along with Witten [11], they proved these on-shell recursion relations using 
just a knowledge of the factorisation properties of the amplitudes and Cauchy's theorem. 
The generality of the proof has led to their application in many diverse areas beyond that 
of massless gluons and fermions in gauge theory [10, 13]. There have been extensions to 
theories with massive scalars and fermions [14, 15, 16] as well as amplitudes in gravity [12]. 

Similarly "on-shell" approaches can also be constructed at loop level. The unitarity 
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of the perturbative ^-matrix can be used to produce compact analytical results by "glu- 
ing" together on-shell tree amplitudes to form the desired loop amplitude. This unitarity 
approach has been developed into a practical technique for the construction of loop ampli- 
tudes [17, 18, 19], initially, for computational reasons, for the construction of amplitudes 
where the loop momentum was kept in D = 4 dimensions. This limited its applicability to 
computations of the "cut-constructible" parts of an amplitude only, i.e. (poly) logarithmic 
containing terms and any associated n 2 constants. Amplitudes consisting of only such terms, 
such as supersymmetric amplitudes, can therefore be completely constructed in this way. 
QCD amplitudes contain in addition rational pieces which cannot be derived using such cuts. 
The "missing" rational parts are constructible directly from the unitarity approach only by 
taking the cut loop momentum to be in D = 4 — 2e dimensions [20] . The greater difficulty 
of such calculations has, with only a few exceptions [21, 22], restricted the application of 
this approach, although recent developments [23, 24, 25] have provided new promise for this 
direction. 

The generality of the foundation of on-shell recursion relation techniques does not limit 
their applicability to tree level processes only. The "missing" rational pieces at one-loop, in 
QCD and other similar theories, can be constructed in an analogous way to (rational) tree 
level amplitudes [26, 27]. The "unitarity on-shell bootstrap" technique combines unitarity 
with on-shell recursion, and provides, in an efficient manner, the complete one-loop ampli- 
tude. This approach has been used to produce various new analytic results for amplitudes 
containing both fixed numbers as well as arbitrary numbers of external legs [28, 29, 30]. 
Other newly developed alternative methods have also proved fruitful for calculating rational 
terms [31, 32, 33, 34]. In combination with the required cut-containing terms [35, 36, 37] 
these new results for the rational loop contributions combine to give the complete analytic 
form for the one-loop QCD six-gluon amplitude. 

The development of efficient techniques for calculating, what were previously difficult 
to derive rational terms, has emphasised the need to optimise the derivation of the cut- 
constructible pieces of the amplitude. One-loop amplitudes can be decomposed entirely in 
terms of a basis of scalar bubble, scalar triangle and scalar box integral functions. Deriving 
cut-constructible terms therefore reduces to the problem of finding the coefficients of these 
basis integrals. For the coefficients of scalar box integrals it was shown in [38] that a 
combination of generalised unitarity [19, 39, 40, 41], quadruple cuts in this case, along with 
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the use of complex momenta could be used, within a purely algebraic approach, to extract 
the desired coefficient from the cut integrand of the associated box topology. 

Extracting triangle and bubble coefficients presents more of a problem. Unlike for the case 
of box coefficients, cutting all the propagators associated with the desired integral topology 
does not uniquely isolate a single integral coefficient. Inside a particular two-particle or triple 
cut lie multiple scalar integral coefficients corresponding to integrals with topologies sharing 
not only the same cuts but also additional propagators. These coefficients must therefore 
be disentangled in some way. There are multiple directions within the literature which have 
been taken to effect this separation. The pioneering work by Bern, Dixon, Dunbar and 
Kosower related unitarity cuts to Feynman diagrams and thence to the scalar integral basis, 
this then allowed for the derivation of many important results [17, 18, 19]. More recently the 
technique of Britto et. al. [23, 24, 25, 35, 36] has for two-particle cuts and the its extension to 
triple cuts by Mastrolia [42] , highlighted the benefits of working in a spinor formalism, where 
the cut integrals can be integrated directly. Important results obtained in this way include 
the most difficult of the cut-constructable pieces for the one-loop amplitude for six gluons 

with the helicity configurations A 6 (-\ 1 1 — ) and A 6 (—-\ h+). The cut-constructible 

parts of Maximum-Helicity- Violating (MHV) one-loop amplitudes were found by joining 
MHV amplitudes together in a similar manner to at tree level [43]. This method has been 
applied by Bedford, Brandhuber, Spence and Travaglini to produce new QCD results [37]. 
In the approach of Ossola, Papadopoulos and Pittau [44, 45] it is possible to avoid the need 
to perform any integration or use any integral reduction techniques. Coefficients are instead 
extracted by solving sets of equations. The solutions of these equations include the desired 
coefficients, along with additional "spurious" terms corresponding to coefficients of terms 
which vanish after integrating over the loop momenta. 

The many-fold different processes and their differing parton contents that will be needed 
at current and future collider experiments suggests that some form of automation, even of 
the more efficient "on-shell" techniques, will be required. From an efficiency standpoint, 
therefore, we would ideally wish to minimise the degree of calculation required for each step 
of any such process. Here we propose a new method for the extraction of scalar integral 
coefficients which aims to meet this goal. The technique follows in the spirit of the simplicity 
of the derivation of scalar box coefficients given in ref. [38]. Desired coefficients can be 
constructed directly using two-particle or triple cuts. The complete one-loop amplitude 
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can then be obtained by summing over all such cuts and adding any box terms and rational 
pieces. Alternatively our technique can be used to extract the bubble and triangle coefficients 
from a one-loop amplitude, generated for example from a Feynman diagram. Hence the 
technique is acting as an efficient way to perform the integration. 

We use unitarity cuts to freeze some of the degrees of freedom of the integral loop mo- 
mentum, whilst leaving others unconstrained. This then isolates a specific single bubble 
or triangle integral topology and hence its coefficient. Within each cut there remain ad- 
ditional coefficients. In the triangle case those of scalar box integrals. In the bubble case 
both scalar box and scalar triangle integrals contribute. Disentangling our desired coefficient 
from these extra contributions is a straightforward two step procedure. First one rewrites 
the loop momentum inside the cut integrand in terms of its unconstrained parameters. In 
the triangle case there is a single parameter, and in the bubble case there are a pair of 
parameters. Examining the behaviour of the integrand as these unconstrained parameters 
approach infinity then allows for a straightforward separation of the desired coefficient from 
any extra contributions. The coefficient of each basis integral function can therefore be 
extracted individually in an efficient manner with no further computation. 

This paper is organised as follows. In section II we outline the notation used throughout 
this paper. In section III we proceed to present the basic structure of a one-loop amplitude 
in terms of a basis of scalar integral functions. We describe in section IV our procedure 
for extracting the coefficients of scalar triangle coefficients through the use of a particular 
loop-momentum parameterisation for the triple cuts along with the properties of the cut 
as the single free integral parameter tends to infinity. Section V extends this formalism 
to include the extraction of scalar bubble coefficients. The two-particle cut used in this 
case contains an additional free parameter and requires an additional step in our procedure. 
Finally in section VI we conclude by providing some applications which act as checks of 
our method. Initially we examine the extraction of various basis integral coefficients from 
some common one-loop integral functions. We then turn our attention to the construction 
of the coefficients of some more phenomenologically interesting processes. These include the 

three-mass triangle coefficient for the six photon amplitude A e ( — I 1 h), as well as 

a representative three- mass triangle coefficient of the process e + e" — > q + qT 9 9 + ■ Finally 
we construct the complete cut-containing part of the amplitude Ag~ loop (l~, 2~, 3 + , 4 + , 5 + ) 
and discuss further comparisons against coefficients of more complicated gluon amplitudes 
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contained in the literature. 
II. NOTATION 

In this section we summarise the notation used in the remainder of the paper. We will 
use the spinor helicity formalism [47, 48], in which the amplitudes are expressed in terms of 
spinor inner-products, 

(jl) = {r\l + )=u-(k j )u + (k l ) ) [jl] = (j + \l-}=u + (k J )u4k l ), (2.1) 

where u±(k) is a massless Weyl spinor with momentum k and positive or negative chirality. 
The notation used here follows the QCD literature, with [ij] = sign(/c°A;°)(j i)* for real 
momenta so that, 

{ij)\ji]=2k i -k j = s ij . (2.2) 
Our convention is that all legs are outgoing. We also define, 

Xi = u + (ki), Xi = u-(ki). (2.3) 
We denote the sums of cyclicly-consecutive external momenta by 

K., = k * +K +1 + --- + kj_ 1 + kj , (2.4) 
where all indices are mod n for an n-gluon amplitude. The invariant mass of this vector is 

= K't, ■ ( 2 - 5 ) 
Special cases include the two- and three-particle invariant masses, which are denoted by 

Sij = K^j = (ki + kj) = 2ki ■ kj, Sijk = (ki + kj + k^ . (2-6) 

We also define spinor strings, 

(i-\(t±t)\j-) = (ta)[aj}±(tb)[bj}, 

+ + = [ia](a-\(f + fi\j-) + [ib](b-\{f + fi)\j-) . (2.7) 
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III. UNITARITY CUTTING TECHNIQUES AND THE ONE-LOOP INTEGRAL 
BASIS 



Our starting point will be the general dimensionally-regularised decomposition of a one- 
loop amplitude into a basis of scalar integral functions [18, 53] 



2\e 



\ i ij ijk / 

The scalar bubble, triangle and box integral functions are denoted by B , C and D respec- 
tively, and along with rr their explicit forms can be found in Appendix C. The bi, Cij and 
dijk are their corresponding rational coefficients. Any e dependence within these coefficients 
has been removed and placed into the rational, TZ n , term. The problem of deriving the 
one-loop amplitude is therefore reduced to that of finding the coefficients of these scalar 
integral functions and any rational terms when working in D = 4 dimensions. 

We are going to consider obtaining these coefficients via the application of various cuts 
within the framework of generalised unitarity [19, 39, 40, 41]. In general our cut momenta 
will be complex, so for our purposes we define a "cut" as the replacement 

^(2tt)5((/ + ^) 2 ). (3.2) 



By systematically constructing all possible unitarity cuts we can reproduce every integral 
coefficient of a particular amplitude. Alternatively, application of the same procedure of 
"cutting" legs can be used to extract from a one-loop integral the corresponding coefficients 
of the standard basis integrals making up that particular integral, in a sense acting as a form 
of specialised integral reduction. This approach follows in a similar vein to that adopted by 
Ossola, Papadopoulos and Pittau [44]. 

The most straightforward implementation of the technique we present here is when the cut 
loop momentum is massless and kept in D = 4 dimensions. Eq. 3.1 therefore contains, within 
the term TZ n , any rational terms missed by performing cuts in only D = 4. Approaches for 
deriving such terms independently of unitarity cuts exist and so we do not concern ourselves 
with these here [23, 24, 26, 27, 29, 30, 31, 32, 33, 34, 44, 45]. 

As was demonstrated in [38], the application of a quadruple cut, as shown in figure 1, to 
^n~ lo ° P uniquely identifies a particular box integral topology D (Kf,Kj,Kl) and hence its 
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FIG. 1: A generic quadruple cut used to isolate the scalar box integral Do(Kf , K? K£ 



coefficient. This coefficient is then given by 

I 2 

dijk = ^ ^ ] ^li}ijk-,a)^2i}ijk\a)^{J'ijk\a)^i{}ijk\a) 
a=l 



(3.3) 



where Ujk- a is the a th solution of the cut loop momentum I that isolates the scalar box 
function D Q (K?, K?, K%), there are 2 such solutions. Eq. 3.3 applies as well to the cases 
when one or more of the four legs of the box is massless. This is a result of the existence, for 
complex momenta, of a well-defined three-point tree amplitude corresponding to any corner 
of a box containing a massless leg. 

Applying a triple cut to the amplitude v4^~ loop does not isolate a single basis integral. 
Instead we have a triangle integral plus a sum of box integrals obtained by "opening" a 
fourth propagator. This can be represented schematically via 



(/' 
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(4vr 



l2-e ^ 



C (Kl K)) + £ d ijk D (Kf, K], K 



(3.4) 



where the additional terms correspond to "opening" the Ki leg or the Kj leg instead of the 
— (Ki + Kj) leg. Similarly in the case of a two-particle cut we again cannot isolate a single 
basis integral by itself. Instead we get additional triangle and box integrals corresponding 
to "opening" third and forth propagators. Schematically this is given by 



rr 



(// 



2\€ 



(47T) 



2-e 



b t B (K?) + «iC {Kf, K)) + ^ d ljk D (K?, i^ 2 ) + . . . ) , (3.5) 

3 jk 



where again the additional terms are boxes with the Ki leg or the Kj legs "opened" . Whilst 
not isolating a single integral each of the above cuts does single out either one scalar triangle, 
in the triple cut case, or one scalar bubble, in the two-particle cut case. Disentangling 
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these single bubble or triangle integral functions from the contributions of the remaining 
basis integrals will allow us to directly read off the corresponding coefficient. Applying all 
possible two-particle, triple and quadruple cuts then enables us to derive the coefficients of 
every basis integral function. 

IV. TRIPLE CUTS AND SCALAR TRIANGLE COEFFICIENTS 

A triple cut contains not only contributions for the corresponding scalar triangle integral, 
but also contributions from scalar box integrals which share the same three cuts as the 
triangle. Of the four propagators of a scalar box integral, three will be given by the three 
cut legs of the triple cut loop integral. The forth propagator will be contained inside the cut 
integrand in a denominator factor of the form (/ — P) 2 , which corresponds to a propagator 
pole. Ideally we want to separate terms containing such poles from the remainder of the 
cut integrand. The remaining term will be the scalar triangle integral multiplied by its 
coefficient for that particular cut. 

The three delta functions of a triple cut constrain the cut loop momentum such that 
only a single free parameter of the integral remains, which we label t. We can express the 
loop momentum in terms of this parameter using the orthogonal null four- vectors, af , with 
% = 1,2,3, specific forms for these basis vectors are presented in section IV A. The loop 
momentum is then given by 



Denominator factors of the cut integrand depending upon the cut loop momentum, can be 
written as propagators of the general form, (/ — P) 2 . When these propagators go on-shell 
they will correspond to poles in t. These poles will be solutions of the following equation 



If we consider t to be a complex parameter then we can use a partial fraction decomposi- 
tion in terms of t to rewrite an arbitrary triple-cut integral. For the extraction of integral 
coefficients we need only work with integrals in D = 4 dimensions. We also drop an overall 
denominator factor of l/(27r) 4 which multiplies all integrals. The partial fraction decompo- 
sition is therefore given, in the case when we have applied a triple cut on the legs I 2 , {l — Ki) 2 



p = a^t + -< + a%. 



(4.1) 



(/-P) 2 = 2{a ■ P)t + 2{ ai ■ P)- + 2{a 2 ■ P) - P 2 = 0. 



(4.2) 
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and (Z - K 2 ) 2 , by 



i=0 

= (2*)*/ dH n <% 2 ) f WM] (0 + E ReS 7_^ 2A3 l , (4-3) 
*=° V poles {J} J / 

where U — l — Ki and Z = I- This is a sum of all possible poles of t, labelled here as the set 
{j}, contained in the cut integrand denoted by AiA 2 A 3 . Pieces of the integrand without a 
pole are contained in the Inf term, originally given in [30] , and defined such that 

lim ([Inf^^a] (t) - A 1 (t) A 2 (t) A 3 (t)) = 0. (4.4) 

t— +oo 

In general [Inf 2 A 1 A 2 A 3 ](t) will be some polynomial in t, 

m 

[Inf t A l A 2 A 3 ](t) = J2te\ (4-5) 

i=0 

where m is the leading degree of large t behaviour and depends upon the specific integrand 
in question. 

After applying the three delta functions constraints we see that taking the residue of 
AiA 2 A 3 at a particular pole, t = t , removes any remaining dependence upon the loop 
momentum. Hence we can write 

/ dH f[ W) ReSt ^ MM ~ lim l(t - t )A 1 A 2 A 3 ] j dH f[ m)T^r- (4-6) 
J i=o J i=0 

Where on the right hand side of this we understand the integral, f dH, as over the parame- 
terised form of / in terms of t and the three other degrees of freedom. In the cut integrand 
the only source of poles in t is from propagator terms of the type 1/(1 — P) 2 . Generally 
each such propagator, when on-shell, contains two poles due to the quadratic nature, in t, 
of eq. (4.2). If we label these solutions t± then we can write a triple-cut scalar box in terms 
of these poles as 

/ d %^jrhf ~i^r{j ^n^'^T " / ^S^'^r) ■ < 4 - 7 > 

From comparing this to eq. (4.6) we see that all residue terms of eq. (4.3) simply correspond 
to pieces of triple-cut scalar box functions multiplied by various coefficients. 
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Therefore we can associate all residue terms with scalar boxes, meaning that our triple 
cut amplitude can be written simply as 

(2vr) 3 J dH f[ SfflA^Az = (2tt) 3 J dtJ t (jr fA + dlD ° Ut - ( 4 - 8 ) 

i=0 ' \ i=0 ' boxes {1} 

This is a sum over the set {/} of possible cut scalar boxes, Dq u1 , and their associated 
coefficients, di, along with a power series in positive powers of t. In eq. (4.8) we have 
integrated over the three delta functions after performing the integral transformation from 
to t, the Jacobian of which, and any additional factors picked up from the integration 
is contained in the factor J t . The limit m of the summation is the maximum power of 
t appearing in the integrand, which in turn is the maximum power of / appearing in the 
numerator of the integrand. In general for renormalisable theories, such as QCD amplitudes, 
m < 3. 

We must now turn our attention to answering the question of what do the remaining 
terms correspond to? To do this we need to understand the behaviour of the integrals 
over positive powers of t. There is a freedom in our choice of the parameterisation of the 
cut-loop momentum. This freedom extends, as we will prove in section IV A, to choosing 
a parametrisation where the integrals over all positive powers of t vanish. Doing this then 
reduces the cut integrand to 

(2vr) 3 J dH f[ 5(l?)A 1 A 2 A 3 = {2-n fh J dtJ t + £ d t D™\ (4.9) 
' " ' boxes 0} 

The remaining integral is now simply that of a triple-cut scalar triangle, multiplied by the co- 
efficient f . For the triple-cut scalar triangle integral, CQ Ut (Kf,Kj), given by — (27r) 3 / dtJ t , 
the triple cut form of eq. (C4), we find that its corresponding coefficient is given simply by 

Qj = -[Inf t A^3]W| t=o , (4-10) 

which is just the first term in the series expansion in t of the cut-integrand at infinity. 

The simplicity of this result relies crucially upon two facts. The first is that on the 
triple cut the integral is sufficiently simple that it can be decomposed into either a triangle 
contribution or a box contribution. This is important as it allows us to easily distinguish 
between the two types of term. As an example consider a linear box which contains a 
numerator factor constructed such that it vanishes at the pole contained in the denominator, 
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but without being proportional to the denominator itself. To which basis integral does this 
term contribute to? In the simplest case such a term would look like 

/ *n«M>© - / "fi'm^S - { ^ I *n««>. 

and hence must contribute entirely to the triangle integral, it contains no box terms. Here we 
have chosen a simplified loop momentum parameterisation in terms of two basis spinors \a + ) 
and \a + ) such that (IP) = t(aP) + (aP). This then contains a pole in t at to — —(aP)/(aP) 
and we have chosen the spinor \W + ) such that (aW) = —to(aW). 

The second crucial fact is the vanishing of the other integrals over t so that the complete 
scalar triangle integral is given by only the remaining integral over t°. Hence the coefficient 
is given by a single term. Furthermore, the use of a complex loop momentum also means 
that we can apply this formalism to the extraction of scalar coefficients corresponding to 
one- and two-mass triangles as well as three-mass triangles. As discussed above for the case 
of box coefficients, this is a result of the possibility of a well-defined three-point vertex when 
using complex momentum, enabling in these cases the construction of non-vanishing cut 
integrands. 

A. The momentum parameterisation 

We wish to compute the coefficient of the scalar triangle singled out by the triple cut 
given in figure 2. The cut integral when written in terms of tree amplitudes is 




FIG. 2: The triple cut used to compute the scalar triangle coefficient of Co (Kf , • 
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r 2 

(2tt) 3 / dH\{ 6(l*)A%™ Ci+2 (-l, Cl , . . . , (c 3 - 1), h)ATc 3+ 2(-h, c 3 , . • • , (c 2 - 1), l 2 ) 

J i=0 

x^ c C2+Cl+2 H 2 ,c 2 ,...,( Cl -l),/), (4.11) 



with lx — I — Ki — I — K Cl C3 _i and l 2 = I — K 2 = I + K C2Cl _i, so that K x = if cl ... C3 _i and 
K 2 = -X C2 ... Cl -i- 

Our first step will be to find a parameterisation of / in terms of the single free integral 
parameter remaining after satisfying all three of the cut delta functions constraints, 

Z 2 = 0, lj = (I - K x f = 0, and l\ = (I - K 2 f = 0. (4.12) 

Each of the three legs can be massive or massless. We will deal with the general case of three 
massive legs explicitly here. The cases with massless legs are then easily found by setting 
the relevant mass in the parameterisation to zero. We will find it very convenient to express 
Z M in terms of a basis of momentum identical to the momenta l\ and l 2 used by Ossola, 
Papadopoulos and Pittau [44]. We will write these momenta in the suggestive notation K\ 
and K\ and define them via 

7 

KY = K% - —K^, (4.13) 

7 

with 7 = (ir{'-|#2|tfi") = (K2~\ft[\ K 2~) and S { = Kf. Each momentum K\, K\ is 
the massless projection of one of the massive legs in the direction of the other masslessly 
projected leg. A more practical definition of K\ and in terms of the external momenta 
alone, can be found by solving the above equations for K\ and K\, so that in terms of Si, 
S 2 , K± and K 2 we have 

M _ K? - {SiHW M _ Kl - (S 2 ft)Kt 

1 " i - (SiS 2 M ' K * - i- (SiS 2 m ■ (4 - 14j 

In addition 7 can be expressed in terms of the external momenta, 

7± = (Ki ■ K 2 ) ± VA, A = (Ki ■ K 2 f - K\Kl. (4.15) 

When using eq. (4.10) we must average over the number of solutions of 7. In the three-mass 
case there are a pair of solutions. For the one- and two-mass cases, when either K\ = or 
K\ = 0, then there is only a single solution. 
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After satisfying the three constraints given by eq. (4.12) we write the spinor components 
of Z M in terms of our basis K\ and K\ as 

(r| = t(Kl>-\+a i(Kl>-\, 

(/+| = ^l(K\' + \ + (K>' + \, (4.16) 



t 



where 



_ Sift-Si) _ ^2(7-^1) ,„ 17 x 

Written as a four-vector, P is given by 

F = a 02 K^ + a 01 KY + liKin-flKln + ^^(K b 2 '-\Y\K b n- (4.18) 
We can also use momentum conservation to write component forms for the other two cut 
momenta k with i = 1,2, 



= t(Kl>-\+ ail (Kl>-\, 

(III = ^(Kl' + \ + (K^\, (4.19) 

where the are given in Appendix A. 

A final point is that after having integrated over the three delta function constraints and 
performed the change of variables to the momentum parameterisation of eq. (4.16) we have 
the factor J t = 1/(27) contained in eq. (4.8). We always associate this factor with the scalar 
triangle integral and so its explicit form does not play a role in our formalism. 



B. Vanishing integrals 

As we have remarked previously, the simplicity of the method outlined here rests crucially 
upon the properties of the momentum parameterisation we have used. The key feature is the 
vanishing of the integrals over t. It can easily be shown that within our chosen momentum 
parameterisation, of section IV A, any integral of a positive or negative power of t vanishes. 
Following an argument very similar to that used by Ossola, Papadopoulos and Pittau [44] 
we use {K^\t,\K^) = 0, (K^\$ 2 \K^) = and (J^Vl^K^W^'*) = 0, to 
show that 

/ d4l {K \li K } nn =0 =► fdtJ t - = forn>l, 
J Plffi J > 

/ dH ^ m^i ^ = ^ / dtJttn = for n ^ L ( 4 - 2 °) 
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The vanishing of these terms then leads directly to our general procedure, encapsulated in 
eq. (4.10), which is to simply express the triple cut of the desired scalar triangle in the 
momentum parameterisation given by eq. (4.16) and then take the t° component of a series 
expansion in t around infinity. 



V. TWO-PARTICLE CUTS AND SCALAR BUBBLE COEFFICIENTS 



In the same spirit as the triangle case we now wish to extract the coefficients of scalar 
bubble terms using, in this case, a two-particle cut. Now a two-particle cut will contain in 
addition to our desired scalar bubble both scalar boxes and triangles, all of which need to 
be disentangled. What we will find, though, is that naively applying the technique as given 
for the scalar triangle coefficients will not give us the complete scalar bubble contribution. 

The reason for this is straightforward to see. A two-particle cut places only two constraints 
on the loop momentum and so we can parameterise it in terms of two free variables, which we 
will label t and y. Consider rewriting the cut integrand in a partial fraction decomposition 
in terms of y. Schematically, therefore, the two-particle cut of the legs I 2 and (/ — Ki) 2 can 
be written as 



{2-nfJ d 4 l fj 5(l 2 ) A, A 2 = {2-KfJ dtdyJ t , y I [Inf^ A 2 ] (y) + £ 

\ poles {j} 



Res y=yj A X A 2 



i=0 



y-Vo 



,(5.1) 



where again {j} is the sum over all possible poles, this time in y, and J t , y contains any 
terms from the change into the parameterisation of y and t as well as any pieces picked up 
by integrating over the two delta functions. So far this seems to be similar to the triangle 
case, but with the residue terms now corresponding to triangles as well as boxes. As we 
have two parameters though we can consider a further partial fraction decomposition, this 
time with t, giving 

(2n) 2 J d 4 lf[ 5(l 2 )A 1 A 2 = 



i=Q 



dtdyj t:y ( [Inf t [Inf,^] (y)} (t) + 
Res t=tl [Ini y A 1 A 2 ] (y) 



+ E 

poles {/} 



t-ti 




+ E 

poles {j},{i} 



, (5.2) 
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where here {/} is the sum over all possible poles in t. The general dependence of the cut 
integral momentum, Z M , on the free integral parameters t and y can be written in terms of 
null four-vectors af with i — 0,1, 2, 3, 4 such that I 2 = 0. An explicit form for these will be 
presented in section VA. We then define / M by 



Again residues of pole terms will correspond to the solutions of (/ — P) 2 = and hence it 
is straightforward to see that the final term of eq. (5.2), containing the sum of residues in 
both y and t, has both of these free parameters fixed. Any such terms must contain at least 
one propagator pole. Also the numerator will be independent of any integration variables, 
as both y and t are fixed. Thus all such terms will correspond to purely scalar triangle and 
scalar box terms. Looking at the second and third terms of eq. (5.2) we might also, at least 
initially, want to associate these terms with contributions to scalar triangle terms only and 
hence naively conclude that only the first term of eq. (5.2) contributes to the scalar bubble 
coefficient. This assumption though would be wrong. 

The crucial difference between the single residue terms of eq. (5.2) and those of eq. (4.3) 
is the parameterisation of the loop momentum which is being used. Taking the residue of 
a pole term at a particular point y freezes y such that we force a particular momentum 
parameterisation upon these triple-cut terms. Importantly, in general this particular forced 
momentum parameterisation is such that the integrals over t in the second and third terms 
of eq. (5.2) now no longer vanish. 

If only scalar triangle contributions came from the integrals over t then this would not 
be an issue; we could just discard these terms as not relevant for the extraction of our 
bubble coefficient. What we find though, through a simple application of Passarino-Veltman 
reduction techniques, is that these integrals contain scalar bubble contributions, B , with 
coefficients b, 



where J' t is the relevant Jacobian for this parameterisation of the loop momentum and c is 
the coefficient corresponding to the scalar triangle contribution, C . We cannot therefore 
simply discard the residue pieces of eq. (5.2), as we could in the triangle case, if we want to 
derive the full scalar bubble coefficient. Furthermore, there is an additional complication. 




(5.3) 




(5.4) 
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We will see that the integrals over powers of y contained in the first term of eq. (5.2) also 
do not vanish in general and hence must also be taken into account. 

There is a limit to the maximum positive powers of y and t that appear in the rewritten 
partial-fractioned decomposition of the integral. For renormalisable theories, such as QCD, 
up to three powers of t appear for triangle coefficients and up to four powers of y for bubble 
coefficients. Therefore the power series in y and t of the Inf operators will always terminate 
at these fixed points. It is then straightforward, as we will discuss in section VD and 
section VB, to derive the general form for all possible non- vanishing contributing integrals, 
over powers of y and t, in terms of their scalar bubble contributions. 

Calculation of the scalar bubble coefficient therefore requires a two stage process. First 
take the Inf^ and lnf t pieces of the cut integrand and replace any integrals over y with 
their known general forms, as we shall see integrals proportional to t will vanish. Secondly 
compute all possible triple cuts that could be generated by applying a third cut to the two- 
particle cut we are considering. To these terms then apply, not the parameterisation we used 
in section IV, but the parameterisation forced upon us by taking the residues of the poles in 
y, which we will derive in section V C. This is equivalent to calculating all the contributions 
from the residues of the partial fraction decomposed cut integrand of eq. (5.2). Within these 
terms we then replace any integrals of powers of t with their known general forms. Finally 
we sum all the contributing pieces together to get the full scalar bubble contribution and 
hence its coefficient. Our final result for assembling the bubble coefficient is then given by 
eq. (5.28). 

A. The momentum parameterisation for the two-particle cut 

We want to extract the scalar bubble coefficient obtainable from the application of the 
two-particle cut given in figure 3. This two-particle cut can be expressed in terms of tree 
amplitudes as 

(2tt) 2 I dH f[ 8{ll)A%°l Ci+2 {-l, ( Cl + 1), . . . , c 2 , ZiK- c C2+Cl+2 (-/i, (ca + 1), . . . , ci, /), (5.5) 

"* i=0 

with Zi = I - K Cl+ i,„ C2 = I - Ki. 

A bubble can be classified entirely in terms of the momentum of one of its two legs, which 
we label K±, and so we will find it useful to express the cut loop momentum / in terms of 
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/ 



FIG. 3: The two-particle cut for computing the scalar bubble coefficient of Bq(K 2 ). 

the pair of massless momenta K\ and x defined via 

K Y = K» - ^V, (5-6) 

7 

here 7 = (x ± \fti\x ± ) = (x^ I 1 X ± ) • The arbitrary vector \ can be chosen independently 
for each bubble coefficient as a result of the independence of the choice of basis representation 
for the cut momentum. In the two-particle cut case we have only two momentum constraints 

I 2 = 0, and l\ = {l- K{) 2 = 0, (5.7) 

and so we have two free parameters which we will label y and t. The loop momentum can 
then be expressed in terms of spinor components as 

(r\ = t {K\'-\ + ^(l-y){ X -\, 
7 

</+| = ^(K^\ + (x + \. (5.8) 



t 



Written as a four-vector is 



P = yK Y + ^(1 - + Uk^IYIx-) + ^7(1 - v)(x-\-f\«tn- (5-9) 
7 2 27 1 

We can also use momentum conservation to write a component form for the other cut 

momentum. We have 

(H\ = (^'-|-ff(x-|, 

= (y-l)(^ + |+t( X + |. (5.10) 

Furthermore after rewriting the integral in this cut-momentum parameterisation and 
integrating over the two delta function constraints we find the following simple result for 
the constant J t , y contained in eq. (5.1), namely J t>y = 1. 
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B. Non-vanishing integrals 



In the case of the scalar triangles of section IV B crucial simplifications occurred as a result 
of our chosen cut momentum parameterisation. Any integral over a power of t vanished, 
leaving only a single contribution corresponding to the desired coefficient. For the scalar 
bubble coefficient things are not quite as simple. 

We can use (i^l^l v±) = as well as (^}' ± | 7 ^|x ± ) (i^'^l y±) = to show that 



dtdy T = 0, 



/<rt<^|p=0 /««*(f)"(l-vr = 0. (5.11) 

Hence the integrals over all positive and negative powers of t vanish, 

Jdtdyt n = forn^O. (5.12) 

Integrals over positive powers of y, contained within the double Inf piece of the first term 
of eq. (5.2), will not vanish. These integrals are straightforwardly derivable with the aid 
of identities involving the four vector = — (Si/t)/^ which satisfies the constraints 
(Ki ■ n) = and n 2 = —Si. It is then possible to show the following relations in D = 4 
dimensions, and remembering that J tjV = 1, 

Hi± ^P = ° - jdtdy^-y)^ 1 ^, 
niL^l = sl m B^ =► J 'dtdy Q-^ =S\ m B^ 

J dH ^jf- = ( 2m + l ) S l mB vv =► ^ / dtdV = (2m + (5 ' 13) 

where B™ v and Bj} v are Passarino-Veltman reduction coefficients, the explicit forms of which 
are not needed. Solving these equations for the integral of y m leads to the result 

/ dtdy y m = [ dtdy for m > 0. (5.14) 

J m+l J 

Contributions to our desired scalar bubble coefficient from the double Inf piece of eq. (5.2) 
therefore come not only from the single constant t°y° term but also from terms proportional 
to integrals of t°y m . This is not the end of the story. As described above, there can be further 
contributions from the second and third residue terms generated in the decomposition of 
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eq. (5.2). We could proceed from the cut integrand to explicitly calculate these residue 
terms. However as we will shall see, a more straightforward approach is to derive these 
terms by relating them to triple cuts. 

C. The momentum parameterisation for triple cut contributions 

We wish to relate the contributions to the bubble coefficient of the residue pieces, sep- 
arated in the decomposition of eq. (5.2), to triple cuts in a specific basis of the cut-loop 
momentum. To find this basis we will apply the additional constraint 

(l + K 2 ) 2 = 0, (5.15) 

to the two-particle cut momentum of section V A. Note that here we label the U K 2 " leg as K 2 
in contrast to (—K 2 ) as we did in the triangle coefficient case of section IV A. This constraint 
corresponds to the application of an additional cut which would appear as 5((l + K 2 ) 2 ) inside 
the integral. This additional constraint, applied to the starting point of the two-particle cut 
loop momentum, forces us to use K\ and \ as the momentum basis vectors of I. Importantly, 
this differs from the basis choice for the triple cut momenta developed in section IV A, which 
leads to the differing behaviour of these triple-cut contributions. 

The presence of y in both and directs us for reasons of efficiency to choose to 
use eq. (5.15) to first constrain y, leaving t free. Looking at eq. (5.9) we see that as is 
quadratic in y then there are two solutions to this constraint, y±, which are given by 

2S 1 {x-\W 2 \K 1 ' ) V J 

±^(s^ X -\t 2 \K\-) + 2t 7 (tfi -K 2 )) 2 -4^27* (*7- (xi^I^Jo)). ( 5 - 16 ) 

On substituting these two solutions into the two-particle cut momentum of eq. (5.8) we 
obtain our desired triple-cut momentum parameterisation. 

Our final step is then to relate the triple-cut integrals defined in this basis to the residue 
terms of eq. (5.2). Rewriting the triple cut integral after the change of momentum parame- 
terisation and integrating over all but the third delta function gives the general form 

(2vr) 3 J dtdyJl(6(y-y + ) + 6(y-y-))M(y,t), (5.17) 
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where A4(y,t) is a general cut integrand and 

A = -, 1 = (5-18) 

\j (s l ( X ~\$2\K'n + 2t 7 (K, ■ K 2 )) - 4S 1 S 2l t (f 7 - |# 2 |tf{'-)) 

Upon examination of a general residue term we find that it corresponds to an integral of the 
form 

{2it) 2 ifdtdyJ t , y Res /?J; y '% = / dtdy J[ (5(y - y + ) + 5(y - y_)) M(y, t){5.19) 

and hence that residue contributions are given, up to a factor of (—1/2), by the triple cut. 

This result applies equally when S 2 = K\ = 0, corresponding to a one or two-mass 
triangle, when the appropriate scale is set to zero in eq. (5.16) and eq. (5.18). The momentum 
parameterisation in this simplified case is contained in Appendix B. 



D. More non- vanishing integrals and bubble coefficients 

There is a direct correspondence between a triple cut contribution and a residue contri- 
bution. The sum of all possible triple cuts, which contain the original two-particle cut, will 
therefore correspond to the sum of all residue terms. We must now examine how such terms 
contribute to the bubble coefficient itself. 

Unlike for the case of triple cut integrands as parameterised in section IV A we will find 
that there are contributions, specifically in this case bubble coefficient contributions, coming 
from the integrals over t. To see this let us investigate the integrals over t in more detail. 
As an example consider extracting the scalar bubble term coming from a two-mass linear 
triangle (with the massless leg K 2 so that S 2 = 0). We would start from a two-particle cut 
which, after decomposing as eq. (5.1), would give 

(2vr) 2 J ' dnf[6^) { 0j^l (5.20) 

= (2vr) 2 | dnf[5^)^- ] =(2nrjdtdyj' t 

The first term of this is clearly not the complete coefficient, and so we need to obtain the 
bubble contribution contained within the second term. Consider reconstructing this term 
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/ Rpe ( y[K\a]+t\xa] 



[K[K 2 \ [KiK 2 \ (f [K\K,\ + \xK 2 \) 



using a triple cut with the cut loop momentum parameterised in a form given by setting y 
equal to its value at the residue of the pole of this second term. This triple cut term is given 
by 



(2tt)H 



d A l\[5{ll){K-\t\a 



i=0 



(27r) 3 i 



/ (<*Tl*il*r>'+f (x-mKtn) , (5.21) 



where we have used the parameterisation of / M given by eq. (B9) and added an extra overall 
factor of i which would come from the additional tree amplitude in a triple cut. 

Of this triple cut integrand only the first, t dependent, term can give anything other than 
a scalar triangle contribution. To derive the result of this integral over t we will, as we have 
done previously, use our parameterisation of the cut momentum, eq. (5.9), to pick out the 
integral as follows 

Using Passarino-Veltman reduction on the single tensor integral on the left hand side of 
this as well as dropping anything but the contributing bubble integrals of our particular cut 
leaves us with the result 



/ 



' (27r) 3 T^'lifil^") 2 



where -Bo ut (-^i) i s the cut form of the scalar bubble integral of eq. (CI). This non-vanishing 
result for the integral over t, in contrast to that of section IV A, is a direct consequence of 
the cut momentum parameterisation forced upon us when taking the residues contained in 
the two-particle cut integrand with which we started. 

On substituting the result of eq. (5.23) into eq. (5.21) we find that we can write eq. (5.20), 
using the bubble integral given in eq. (CI), as 



lK7\TtAn-\ 

B o ( K i)i (5-24) 



(K2 \$i \K 2I 

which is the known coefficient of the scalar bubble contained inside the linear triangle. 
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Of course, if we had chosen x — K 2 from the beginning, then the first term on the left 
hand side of eq. (5.20) would have been the complete bubble coefficient. In general, if we 
are able to rewrite a two-particle cut integrand such that each term contains only a single 
propagator then we can always choose a different \ — K\, defined via 



KY = K% — ^ <—K^, (5.25) 



for each term individually such that there are no contributions from the residue terms. 
Whether this is both feasible and a more computationally effective approach than calculating 
the residue contributions through the use of triple cuts would depend upon the cut integrand 
in question. 

In general we will be considering processes which contain terms with powers of up to t 3 , 
so we will need to know these integrals. Again these can be found using a straightforward 
application of tensor reduction techniques. When all three legs in the cut are massive these 
integrals over t are given, after dropping an overall factor of l/(27r) 3 witch always cancels 
out of the final coefficient, by 



T(j)=Jdtj t * = {-) Xj r gc,,—— 



B™(Kl). (5.26) 



J=i 

Simply taking the relevant mass to zero gives the forms in the one and two mass cases. A 
was previously defined in eq. (4.15) and we have 

Cii - 

c - 3 c - 3 

(-21 — — g, L-22 — — -, 

Also for later use we define T(0) = 0. 



E. The bubble coefficient 



We have now assembled all the pieces necessary to compute our desired scalar bubble 
coefficient, bj, corresponding to the cut scalar bubble integral BQ Ut (Kj). It is given in general 
not as the coefficient of a single term but by summing together the t°y m terms from both the 
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double Inf in y followed by t as well as residue contributions which we derive by considering 
all possible triple cuts contained in the two-particle cut. The coefficient is given by 



where T(j) is defined in eq. (5.26) and the sum over the set {C tvl } is a sum over all triple 
cuts obtainable by cutting one more leg of the two-particle cut integrand A\A 2 . 

When computing with eq. (5.28) there is a freedom in the choice of x- A suitable choice of 
which can simplify the degree of computation involved in extracting a particular coefficient. 
Particular choices of x can eliminate the need to calculate the second term of eq. (5.28) 
completely, as discussed in section VD. We also note that there are choices of x which 
eliminate the need to evaluate the first term of eq. (5.28), so that the coefficient comes 
entirely from the second term of eq. (5.28) instead. 

VI. APPLICATIONS 

To demonstrate our method we now present the recalculation of some representative 
triangle and bubble integral coefficients. We also discuss checks we have made against other 
various state-of-the-art cut-constructable coefficients contained in the literature. 

A. Extracting coefficients 

To highlight the application of our procedure to the extraction of basis integral coefficients 
we consider deriving the coefficients of some simple integral functions which commonly 
appear, for example, in one-loop Feynman diagrams. 

1. The triangle coefficient of a linear two-mass triangle 

First we consider deriving the scalar triangle coefficient of a linear two-mass triangle with 
massive leg Ki, massless leg K 2 , and a and b arbitrary massless four- vectors not equal to 
K 2 . This is given by the integral 




i [Inf t \Ln£ v A 1 A 2 ] (y)} (t) 



t->0, y m ~^^rr 

' * m+l 



\ £ [Inf t AA 2 A 3 ] (t) 



(5.28) 



{c tri } 




(6.1) 
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Extracting the triangle coefficient requires cutting all three propagators of the integrand. 
We do this here by simply removing the "cut" propagator as we are interested only in the 
integrand. This leaves only 

(a"|/|6->. (6.2) 

Rewriting this integrand in terms of the parameterisation of eq. (4.16) gives 

{a 01 (a-\t 2 \b-)+t(aK\)[ X b]). (6.3) 

As S 2 = we see that ctoi = Si/j and that 7 = 2{K\ • K 2 ). Then taking the t° component 
of the [Inf t ] of this in accordance with eq. (4.10) leaves us with our desired coefficient 

which matches the expected result. 

2. The bubble contributions of a three-mass linear triangle 

Consider a linear triangle with in this case three massive legs, so now K 2 is massive but 
again a and b are arbitrary massless four-vectors, 

J l 2 {l - Ki) 2 (l + K 2 ) 2 ' (6 - 5) 
Extracting the bubble coefficient of the integral B (Kf) is done by cutting the two propa- 
gators I 2 and {I — Ki) 2 . Again cutting the legs is done by removing the relevant propagators 
from the integrand so that it is given by 

£m. (6.6) 

As this contains a single propagator, and therefore a single pole, we could choose to set 
X = K\ (as defined in eq. (5.25)), before performing the series expansions in y and t. For 
this choice of x the bubble coefficient comes entirely from the two-particle cut. Using the 
first term of eq. (5.28) gives directly 



7 2 - S,S 2 7 2 - SiS 
where 7 = (K^'lI^ilK^ ~), a result which is equivalent to the expected answer. 
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In order to demonstrate the procedure of using triple cut contributions in extracting a 
bubble coefficient we will now reproduce this by assuming x 7^ -^2- F° r this case t ne first 
term of eq. (5.28) then gives 

which upon choosing x — a vanishes and so the complete contribution will come from the 
triple cut pieces of eq. (6.5). Cutting the remaining propagator in eq. (6.6) gives us the 
single triple cut term which will contribute. The integrand of this is given, after multiplying 
by an additional factor of % which would come from the third tree amplitude if this was a 
triple cut, by 

% ({al)[lb} + (al)[lb] ) = i {(y + + y.)(aK\)[K\b\ + 2t{aK\)[ab]) , (6.9) 
V y=y+ y=y-J 

where we have set x — a - From eq. (5.16) we have 



7 ((X}'-|# 2 |X}'-) - Si<a-|# 2 |a-)) * + S 1 (a-\^ 2 \K['-) 



y+ + y- = — — ITfr ,^_; , (6.10) 



Hence taking the [Inf t ] of the cut integrand, eq. (6.9), and dropping any terms not propor- 
tional to t leaves 

«(q-|*,IO ( «ti-\*,\K*-) - S.ia-mia-} + 2 M ) (6 n) 

V S,(a-\$ 2 \K\~) \K\b\)' 

which after inserting the result for the t integral given by eq. (5.26) and substituting this 
into the second term of eq. (5.28) gives for our desired coefficient 

.(a-^xlft-) / /r u |r> . b . S 1 . 2S X . L [ah] 



—i- 



4A 



^(K\n^\Kln-^(a~ma~) + ^(a-|^ 2 |^'-> j^) 



= -i^L ((^ • K^aT^b-) - S 1 (a-\^ 2 \b-)) , (6.12) 
where A was given in eq. (4.15). This matches both the expected result and eq. (6.7). 



B. Constructing the one- loop six-photon amplitude Aq(1 ,2 + ,3 ,4 + ,5 ,6 + ) 

Recently an analytic form for the last unknown six-photon one-loop amplitude was ob- 
tained by Binoth, Heinrich, Gehrmann and Mastrolia in ref. [46]. This result was used 
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to confirm a previous numerical result [50]. More recently still further corroboration has 
been provided by [45]. Here we reproduce, as an example, the calculation of the three-mass 
triangle and bubble coefficients, again confirming part of these results. 

Firstly it is a very simple exercise to demonstrate by explicit computation that all bubble 
coefficients vanish. If we were to use the basis of finite box integrals, as defined in [35] , then 
there is only a single unique three-mass triangle coefficient, a complete explicit derivation 
of which we now present. Starting from the cut in the 12 : 34 : 56 channel shown in figure 4 
we can write the cut integrand as 

i 5 

■ K 

s 

/ 

/ 

3" A:! 

FIG. 4: Triple cut six-photon amplitude in the 12 : 34 : 56 channel. 




1QA 4 (~1;\ 1", 2 + , 1%)A A {-1^\ 3", 4 + , l^)A 4 (-l^\B', 6 + , Zj), (6.13) 

with all unlabelled legs photons and l\ = I — and l 2 = I + K 12 . The overall factor of 
16 comes from the differing normalisation conventions between QCD colour-ordered ampli- 
tudes and QED photon amplitudes. Both helicity choices h — hi — h 2 — ± give identical 
contributions. Written explicitly, eq. (6.13) is 

oo, (/l) 2 (/ 2 3) 2 (/i5) 2 

(Z2)(Z 2 2)(Z 1 4)(Z 2 4)(Z6)(Z 1 6)' 1 ' ) 

After inserting the momentum parameterisation of eq. (4.16) this becomes 

1 

32z 

(t{K\2) + a i(A^2» (t{K\2) + a 2l {K\2)) (t{K\A) + a u (^4)) 

(t(K>l) + aoi(^l)) 2 jt(K\3) + a 21 (^3)) 2 (t(K{5) + a n (K b 2 5)) 2 
(t(K\4) + a 2l (K^)) (t{K{6) + a 01 (K*6)) (t(K[6) + a u (K b 2 6))' 

Applying eq. (4.10) implies taking only the t° piece of the [Inf 4 ] of this expression. Averaging 

over both solutions leaves us with our form for the three mass triangle coefficient 

_ m Y (K\imism^ 
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where K\ depends upon the form of y± as given in eq. (4.14). Numerical comparison with 
the analytic result of [46] shows complete agreement. 



C. Contributions to the one-loop Ag(l+,2^.,3 , 4 + ;5 F ,6+) amplitude 

This particular amplitude was originally obtained by Bern, Dixon and Kosower in [19]. 
Making up this amplitude are many box, triangle and bubble integrals along with rational 
terms. Here we will recompute one particular representative three-mass triangle coefficient 
in order to highlight the application of our technique to a phenomenologically interesting 
process. 

Following the notation of [19], we wish to calculate the three-mass triangle coefficient of 
If m (si4, s 23 , s 56 ) = C (si4, s 56 ) of the F cc term. The only contributing cut is shown in figure 
5. We begin by writing down the triple cut integrand for this case 




FIG. 5: Triple cut in the 14 : 23 : 56 channel. 



M-ly?, , fe)^ 4 (-/ 2 ^ 2 , 4 + , 1+ l h g )A,(-l-\ 2r, 3", J*), (6.17) 

where h = I — K 2 z and h = I + K u . Only when h = — , hi = + and h 2 = + do we get a 
contribution. It can be written explicitly as 

(/ 2 5) 2 (// 2 ) 2 (23) 2 



(14)(56)(4/ 2 )(2/)(// 1 )(M 2 )' 
Rewriting this in terms of the loop momentum parametrisation of eq. (4.16) gives 

7 (^i b 5) + « 21 (^5)) 2 (23) 2 

% 

S23 (l - "f) (14) (56) (t(4^> + a 21 (AKD) (t{2K\) + a 01 <22^» ' 



(6.18) 



(6.19) 
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The two solutions of 7 are given by 7± = — (if 23 ■ if 14) ± a/ (if 23 • ifj.4) 2 — S23S14, the a^'s 
are given in Appendix A. 

The application of eq. (4.10) involves taking [Inf t ] of eq. (6.19), dropping all but the t° 
component of the result and then averaging over both solutions of 7 giving the coefficient 

-i E ? T' m • (6.20) 

2 ^S 1 (l-f)(14)(56)(4A-!)(2A'f) 

where again if \ depends upon 7±. Numerical comparison against the solution for this 
coefficient presented in [19], 

l/Vo-l if it . _|t;+\2 _ /ot;\2 c . . c _„^ 

+ ffip, (6.21) 

shows complete agreement, where the operation flip is defined as the exchanges 1 <-> 2, 
3 <-> 4, 5 <-> 6, (ab) <-> [a6]. 

The remaining triangle and bubble coefficients can be derived in an analogous way. We 
have computed a selection of these coefficients for A 6 (l+, 2r ? 3~, 4 + ; 5^, 6+), along with 
coefficients of other amplitudes given in [19], and find complete agreement. 



i [14] «2-|# 14 #23| 


5+) 2 - (25) 2 Sl4S23 ) 


2 (14)[23](56)(2-|if 14 : 


3->(2-|^3 4 | 


1-) 



D. Bubble coefficients of the one-loop 5-gluon QCD amplitude ^5(1 ,2 ,3 + ,4 + ,5 + ) 

This result for the 1-loop 5 gluon QCD amplitude A 5 (l~, 2", 3 + , 4 + , 5 + ) was originally 
calculated by Bern, Dixon, Dunbar and Kosower in [18]. It contains neither box nor triangle 
integrals, only bubbles. We need therefore only compute bubble coefficients. There are only 
a pair of such coefficients, with masses S23 and S234 = ssi- 

For the first cut in the channel if 1 = if 23 we have, for the sum of the two possible helicity 
configurations, the two-particle cut integrand 

2 (l/ 1 ) 2 (l/)(2/)(2/ 1 ) 2 



(23) (45) (51) (4/^(3/!) (//x) 2 ' 
and for the second, in the channel ifi = if 234, 

2 (l/ 1 ) 2 (l/)(2/)(2/ 1 ) 2 



(6.22) 



(6.23) 



(23) (34) (51) (Ah) (5k) (lh) 2 ' 

Focus upon the ifi = if 23 cut initially. There are two pole- containing terms in the 
denominator of this cut. We could choose to partial fraction these terms and then pick x — 
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K 2 in each case to extract the coefficient. Instead though we will derive the coefficient using 
triple cut contributions. Choosing \ = k\ so that after inserting the cut loop momentum 
parameterisation of eq. (5.8) the cut integrand becomes 



2 7 2 (liT 



I) (<2^> - f f (21)) (t(2K\) + f (1 - y) (21) 



^ 2 (23) (45) (51) (jpjfy _ (31)) (<4tf{) - f f (41)) 



(6.24) 



and hence produces no [Inf^flnf^] term. Consequentially the two-particle cut contribution 
to the bubble coefficient vanishes. The same choice of \ similarly removes all two-particle 
cut contributions in the channel K x = K 2U from the corresponding scalar bubble coefficient. 

Examining the triple cuts of the bubble in the K 23 channel shows only two possible 
contributions, again after summing over both contributing helicities, given by 

2i [3/][3/ 2 ](l/ 1 )(l/) 2 (2/ 1 )(2/) 



(45) (51) {IhHhkWkW*) 



(6.25) 



when K 2 = k 3 and 



2i [4/][4/ 2 ](l/ 1 )(l/ 2 ) 2 (2/ 1 )(2/) 2 



(23) (51) (//i)(/i/ 2 )[// 2 ](5/ 2 )(3/) ' 

when K 2 = k^. In both cases K 2 is massless and is of positive helicity so we use the 
parameterisation of the triple cut momenta for y + given in eq. (B2). Then along with 
setting x = ki gives for the first triple cut integrand 



(13)(34)(45)(51)Sf V (IK\) S, (13) J\ (13) 1 " 7 
and for the second 

2*(1^) 2 (24) 2 7 2 /7t (4-1^14-) (llj^f^V W (24) , giMV 6 . 28) 



(23) (34) (45) (51) (14) S 2 \S, (14) (iff*) A(14) N ' 7 

Applying these integrands to the second term of eq. (5.28) by taking [Inf t ], dropping any 
terms not proportional to t and then performing the substitution t % — > T(i) gives for the 

1 ^trc 
3 A 5 



coefficient of the first triple cut simply \A^ CC , and for the second triple cut 



^tree ^+| 


W 23 


|1 + ) 2 


s? 2 (4-|^ 23 |- 


I-) 2 



2(1+] 


;« 2 3 


Jl + > 


3 (4-|^23|- 


4-> 



S 12 " o / 7Z 7 J T7Z\ • ( 6 - 29 ) 



After following the same series of steps as above for the second bubble coefficient with 
Ki = K 2U we find only a single triple cut contributing term corresponding to K 2 = fc 4 . This 
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is related to the second triple cut coefficient derived above via the replacement K 2 3 — > -^234 
and swapping the overall sign. 

After combining the three triple cut pieces above we arrive at the following form for the 
cut constructable pieces of this amplitude 



-A 5 B (s 2 3) 



(4vr) 2 - e V3 



A trec ^ + 


M#23 


|1 + ) 2 


42 (4-|^23|- 


4-) 2 



s u~l ^ytf 2 ^ ) (Bo(s2 3 )-B (s 2U )) ), (6.30) 



3 (4-|#23|4- 

which can easily be shown to match the result given in [18]. 

While this example is particularly simple we have also performed additional compar- 
isons against other results in the literature. Such tests include the cut constructible pieces 
of all two-minus gluon amplitudes with up to seven external legs, originally obtained in 
[18, 37]. Additionally we find agreement for the case when, with six gluon legs, three 
are of negative helicity and adjacent to each other and the remainder are positive helic- 
ity, which was originally obtained in [49]. We have also successfully reproduced the known 
three mass triangle coefficients in Af = 1 supersymmetry for A e (l~, 2 + , 3~, 4 + , 5~, 6 + ) and 
A e (l~, 2~, 3 + , 4~, 5 + , 6 + ), originally obtained in [35]. 



VII. CONCLUSIONS 



The calculation of Standard Model background processes at the LHC requires efficient 
techniques for the production of amplitudes. The large numbers of processes involved along 
with their differing partonic makeups suggests that as much automation as possible is de- 
sired. In this paper we have presented a new formalism which directs us towards this goal. 
Coefficients of the basis scalar integrals making up a one-loop amplitude are constructed in 
a straightforward manner involving only a simple change of variables and a series expansion, 
thus avoiding the need to perform any integration or calculate any extraneous intermediate 
quantities. The main results of this paper can be encapsulated simply by eq. (4.10) and 
eq. (5.28) along with the cut loop momentum given by eq. (4.16), eq. (5.8) and eq. (5.16). 

Although this technique has been presented mainly in the context of using generalised 
unitarity [19, 39, 40, 41] to construct coefficients, and hence the cut-constructible part of 
the amplitude, it can also be used as an efficient method of performing one-loop integration. 
Using the idea of "cutting" two, three or four of the propagators inside an integral, we 
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isolate and then extract scalar basis coefficients. This procedure then allows us to rewrite 
the integral in terms of the scalar one-loop basis integrals, hence giving us a result for the 
integral. 

Different unitarity cuts isolate particular basis integrals. For the extraction of triangle 
integral coefficients this means triple cuts and for bubble coefficients we use a combination 
of two-particle and triple cuts. Extracting the desired coefficients from these cut integrands 
is then a two step process. The first step is to rewrite the cut loop momentum in terms of a 
parameterisation which depends upon the remaining free parameters of the integral after all 
the cut delta functions have been applied. Triangle coefficients are then found by taking the 
terms independent of the sole free integral parameter as this parameter is taken to infinity. 

Bubble coefficients are calculated in a similar if slightly more complicated way. The pres- 
ence of a second free parameter in the bubble case means that we must take into account, not 
only the constant term in the expansion of the cut integrand as the free integral parameters 
are taken to infinity, but also powers of one of these parameters. The limit on the maximum 
power of Z M appearing in the cut integral restricts the appearance of such terms and hence 
we need consider only finite numbers of powers of these free parameters. Additionally it can 
also be necessary to take into account contributions from terms generated by applying an 
additional cut to the bubble integral. The flexibility in our choice of the cut-loop momentum 
parameterisation allows us to directly control whether we need compute any of these triple 
cut terms. Furthermore we can control which of these triple cut terms appears, in cases 
when their computation is necessary. 

As we consider the application of this procedure to more diverse processes than those 
detailed here, we should also investigate the "complexity" of the generated coefficients. 
In the applications we have presented we can see that we produce "compact" forms with 
minimal amounts of simplification required. This is important if we are to consider further 
automation. The straightforward nature of this technique combined with the minimal need 
for simplification means that efficient computer implementations can easily be produced. As 
a test of this assertion we have implemented the formalism within a Mathematica program 
which has been used to perform checks against state-of-the-art results contained in the 
literature. Such checks have included various helicity configurations of up to seven external 

gluons as well as the bubble and three-mass triangle coefficients of the six photon A 6 ( — -\ h 

— h) amplitude. In addition representative coefficients of processes of the type e + e~ — > qqgg 
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have been successfully obtained. 

Our procedure as presented has mainly been in the context of massless theories. Funda- 
mentally there is no restriction to the application of this to theories also involving massive 
fields circulating in the loop. Extensions to include masses should require only a suitable 
momentum parameterisation for the cut loop momentum; the procedure is then expected to 
apply as before. 

In conclusion therefore we believe that the technique presented here shows great potential 
for easing the calculation of needed one-loop integrals for current and future colliders. 
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APPENDIX A: THE TRIPLE CUT PARAMETERISATION 

In this appendix we give the complete detail of the triple cut parameterisation along with 
some other useful results. The three cut momenta are given by 

(l;\=t(K b r\+a tl (K^\, = ^f(K^ + \ + (K^ + \, (Al) 



with 



Si (7 - s 2 ) S 2 (7 - Si t 

a 0l — T^) n n \ 1 ft 02 



(7 2 -SiS 2 )' VZ (7 2 -SiS 2 )' 

Si S 1 S 2 (1 - (fr/7)) _ 7(^2 ~ 7) 

ocu — aoi — 5 777; , 0-12 — ao2 — J- — —5 r, r, , 

7 Y — bik> 2 T - ^1^2 

«2i = a 01 -l = — — — , a 22 = a 02 = = ^ A A2 J 

7 - ^1^2 7 7 ^W2 

along with the identities aoi«o2 = Q;ii«i2 and a ia 02 = a 2i a 22 . When written as four-vectors 

the cut momentum are given by 

It = a l2 KY + oak? + t-iK^iriKi"-) + ^-{^'W^i')- (A3) 
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From these parameterised forms we have the following spinor product identities 

= - y<^>, 
jt 



[Ik] 


OL\1 - 
t 




(111) 


= t(a u 


-a Ql ){K\Kl) 


[M 


a 22 - 
t 


% M = 


(ik) 
ihk] 


= t(a 2 i - 

t 


-a 0l )(K\K\) 
^[K\K\\- 



t \ 7 



(hk) = f(«n - a 21 )(K\Kl) = -t ( 1 - ^ ) (fffffj). (A4) 



and we note that 



-(i-f) (l-^)7 = -7-^ + 5i + ^ = (^i-^) 2 = 53, (A5) 
and so with / = l we have = Si+j, as expected. 



APPENDIX B: THE TRIPLE CUT BUBBLE CONTRIBUTION MOMENTUM 
PARAMETERISATION WHEN K\ = 

In this appendix we give the forms for the triple cut momentum of section VC in the 
case when £2 = 0, i.e. we have a one or two mass triangle. Firstly in these cases the K 2 leg 
is attached to a three-point vertex and so the amplitude for this will contain either [K 2 l] or 
(K 2 l) depending upon the helicity of K 2 . This means that in the positive helicity case only 
the delta function solution 5(y — y + ) survives and for a negative helicity K 2 the 5(y — y_) 
survives. We have for both solutions 

A = 7 r- (Bl) 

(S^x-W^-) + tl(K 2 -\tAK;)) 

The momentum parameterisation for the y + solution is given in spinor components by 

«1 - gg<A- 2 -|, <n^'V|^™, (B2, 
and as a 4-vector by 

" = IS> (dk> <A ' 2ly * |A ' 2+> + ™ Ki ^) ■ (B3) 
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The other momenta are given by 

(ill = t { r^\( K 2\ - -(x~l (it\ = -^r (k;\$i, 

(/2 1 " W (i " 2 {k 1 = "T^T " si W 2 (B4) 

where we have moved the overall factor of t from (ir| to to avoid the presence of a 1/t 
term for aesthetical reasons. The spinor products formed from these are given by 

(«i> = ~(XK\), [Ik] = [K{ X i 

7 

(«2> = o, [ih\ = -j^mi 

(XK\) 

(hl2) = --(xKi), [hh] = -^[K\ X ]. (B5) 

7 &l 

and we see that again, as expected, with / = l , we have (hlj)[ljli\ = Si + j. As we have 
massless legs some spinor products will consequentially vanish. In the two-mass case these 
are 

(ZZ 2 ) = 0, {IK 2 )=0, [l 2 K 2 ] = 0, (B6) 

and for the one-mass case 

[ZiZ 2 ] = 0, (ll 2 ) = 0, (IK 2 ) = 0, (hK 2 } = 0, [hK 3 ] = 0, [l 2 K 3 ] = 0, (B7) 

where K 3 is the momentum of the third leg. 

The momentum parameterisation for the y_ solution is given in spinor components by 



^WM {KMl+s i {x ' 1 <i+l = & {Ktl m 



and as a 4-vector by 



- S (wm^ 1 *^ + 7 <xlyw> ) ■ (B9 » 

The other momenta are given by 

KI = l^^ + l*3 + p4|W!^, «l = g§<*?|. (Bio) 
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The spinor products formed from these are given by 

(Mi) = -ixKi), [Ik] = [K\xl 

bcK\\ 

(hl 2 ) = ^(Klx), [kk] = [XK\]. 

7 



(Bll) 



and again (klj)[ljh] = S i+ j as expected. The vanishing spinor products in the two mass case 
are 



[ll 2 ] = 0, [IK 2 ] = 0, [l 2 K 2 ] = 0, 



(B12) 



and in the one mass case 

(hk) = 0, [lk] = 0, [IK 2 ] = 0, [l 2 K 2 ] = 0, (hK 3 ) = 0, (l 2 K 3 ) = 0. (B13) 



APPENDIX C: THE SCALAR INTEGRAL FUNCTIONS 



The scalar bubble integral with massive leg K\ given in figure 6 is defined as 




FIG. 6: The scalar bubble integral with a leg of mass K\. 



Bo(Kf) = (-»)(4tt 



,2-e 



d 4 ~ 2 n i 



and is given by 
with 



r( | , v) (-^)- e = rr [~ ~ H-Kf) + 2) Oir ) 



r(i + e )r 2 (i-e) 



(CI) 



(C2) 



(C3) 



r(l-2e) 

The general form of the scalar triangle integral with the masses of its legs labelled Kf, 
K\ and K\ given in figure 7 is defined as 
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FIG. 7: The scalar triangle with its three legs of mass Kf, K\ and K^. 



C (Kf,K 2 2 )=t(A7r 



,2-e 



d A ~ 2 n 



1 



(C4) 



(2tt) 4 ^ P(l - Kif(l - K 2 f 

and separates into three cases depending upon the masses of these external legs. In the one 
mass case we have K\ = and K\ = and the corresponding integral is given by 

1 H-Kl) \n 2 (-Kj) 



rr 



e* ■ (-K() \e 2 e ' 2 

If two legs are massive the integral, assuming K\ = 0, is given by 
r r {-Kl)-*- 



+ 0(e), (C5) 



Kl)-* 



C (Kl Kl, - - { _ Kl) _ { _ Kl) 



In (-Kl) - In (-Kl) In 2 (-Kl) - In 2 (-Kl 



+ 



(-Kl)-(-Kl)\ e ' 2 

Finally if all three legs are massive then the integral is as given in [53, 54] 

3 



Co(KlKl) 



Li 2 



1 + iSj 
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Li 2 



1 + iSi 



where 
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_ K\ - K 2 2 - (Ki + K 2 ) 
-Kl + Kl - (Ki + K 2 ) 
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+ 0(e), (C7) 



(C8) 
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with A given by eq. (4.15). 

The general form for a scalar box function is given by 

d„(ki ki KD = H) ( 4,r-< / { C^P(i-K 1 n-K 1 n-K,r (C10) 

The solution of this integral is split up into classes depending upon the masses of the external 
legs. These solutions are labelled as zero mass J" 7 ™, one mass I\ m , two mass hard, I% mh , two 
mass easy 7f me , three mass 7| m and four mass If 71 integrals. The results for which can be 
found in the literature, for example in [17]. 
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